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The development of energy harvesting systems based on fluid/structure interactions is part of 
the global search for innovative tools to produce renewable energy. In this paper, the possibility 
to harvest energy from a flow using vortex-induced vibrations (VIV) of a tensioned flexible cable 
is analyzed. The fluid loading on the vibrating solid and resulting dynamics are computed using 
- - - an appropriate wake-oscillator model, allowing one to perform a systematic parametric study of 

CN ^ the efficiency. The generic case of an elastically-mounted rigid cylinder is first investigated, before 

considering an infinite cable with two different types of energy harvesting : a uniformly spanwise 
distributed harvesting and then a periodic distribution of discrete harvesting devices. The maximum 
harvesting efficiency is of the same order for each configuration and is always reached when the solid 



o 

(N 



^ . body and its wake are in a frequency lock-in state. 



I. INTRODUCTION 

J>^' Geophysical flows such as wind, oceanic currents or tides represent a clean source of energy that is widely available. 

' O , Many devices exist, or are being imagined, to harvest part of this energy, ranging from wind farms to wave energy 

I— { ■ converters or marine turbines, to cite a few. 

j^^ , Among all these, a specific class is based on flow-induced vibrations of structures. Such vibrations may be caused 

• ■ by several distinct mechanisms, see for instance reviews in Blevins [Ij or Naudascher and Rockwell Qj. Fluid elastic 

fj [ instabilities, a major mechanism, have been considered in the perspective of energy harvesting by Tang et al. [3|, 

•'H - Barrero-Gil et al. [J], Doare and Michelin [5] or Singh et al. fis], on plates, cylinders or squares. Vortex- induced 

^.,^' vibrations (VIV) is another important cause of flow-induced vibrations and originates in a strong dynamical coupling 

r^ , between an oscillating body and its fluctuating wake, causing lock-in of frequencies and high-amplitude motion [7| . 

rS i' Bernitsas et al. [8] proposed an energy harvesting device using VIV and the idea has since then been further developed, 

' — ': UM- 

I ■ A key question for all energy harvesting devices is how to access large quantities of energy : because of the low 

^ ' energy density in geophysical flows, large systems are required. In the specific domain of VIV, such large systems 

\^Q i have actually been extensively studied for offshore engineering issues such as mooring cables or risers, see for instance 



C^ Baarholm et al. ll|. This opens new perspectives in harvesting energy, using VIV of very long tensioned flexible 
CN , structures instead of elastically supported rigid short cylinders. 

^— ^ ' Although the local interaction mechanisms between each body section and its wake are still the same, a new 

CO dimension is introduced when considering the spanwise behaviour of such slender structures. This results in a much 
^D ' more complex dynamics, involving vortex-induced waves, travelling or stationary. These have been observed on a 
vN I variety of systems [12h15| | . Extensive numerical computations |16l4l9l | have shown rich dynamics for this configuration 
of a cable coupled with its wake. A simpler approach, using a wake oscillator model, recently allowed to derive some 
generic results on complex issues such as mode switching, time sharing and space sharing in complex flows 
references within]. This approach has been extensively compared with experiments and DNS computations [13 . 
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^ and used for predictive analysis [22|, |23,j . 

j^ , In the present paper we seek to establish some fundamental results on generic conflgurations of energy harvesting 

from VIV of tensioned cables. To this end, the wake oscillator model introduced in Facchinetti et al. [2J] is used, as 
well as its linearized version [2^ . In Section [III the elementary 2D case of an elastically supported rigid cylinder is 
considered. The infinite tensioned cable with continuously distributed harvesters is solved in Section Hill Finally in 
Section HVl the influence of the harvesters distribution is addressed. The general features of this new idea of harvesting 
energy from vortex-induced waves in cables are then discussed in Section [V] 
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II. THE ELASTICALLY MOUNTED RIGID CYLINDER 
A. Fluid-solid model 
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FIG. 1: Energy harvesting from the vortex-induced vibrations of (a) the 2D generic case of an elastically mounted rigid cylinder, 
(b) an infinite tensioned cable with continuous energy harvesters, (c) an infinite tensioned cable with localized energy harvesters. 



We consider first the classical generic two dimensional case (2D) of a cylinder mounted on an elastic and damped 
support, Figure[lja). Let r and h be the damping and stiffness coefficients of the support, D, mg and Y the diameter, 
mass per unit length and cross-flow displacement of the cylinder and p, U the fluid density and velocity. Following 
previous authors [1,, ,24] , the cross- flow dynamics of the cylinder under flow is described by 
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rus + ^pD^Caio 






r^\pDUC,Y^ 



hY = K 



wake-! 



where the fluid loading has been split into three parts : the added mass with a coefficient Cmq 
with the drag coefficient Cd and the remaining part denoting the wake force F^ake- 

Two frequencies exist in the system : (i) that of the cylinder in still fluid {U = 0), namely uj_ 



(1) 
the added damping 



^/h/rrit^ where 



TTit — rUs + pttD'^Cmo/4: is the total mass per unit length, and (ii) the vortex shedding frequency behind the still 
cylinder, ujj = 27rStU/D, where St is the Strouhal number |l[. Using the latter, the dimensionless form of ((T|) reads 



y+[( + -]y + s^y = fwake 



(2) 



where y = Y/D, t — ujfT, 6 = c^s/'-^f and ( ) denotes derivation with respect to the dimensionless time t. The 
dimensionless damping is composed of the support damping ^ = r/mtio f and of the fluid damping j/p where 7 — 
Cu/i-KSt is the stall parameter and p = mt/pD^ is the mass ratio |24ll26|. 

From the point of view of the fluid-solid system, energy harvesting induces a loss of energy, which will be represented 
in the remainder of the paper by the support damping ^. The efficiency of the harvesting may then be defined as in 
Barrero-Gil et al. [1,^] or Zhu et al. |23]. 
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(3) 



where ( ) denotes averaging in time. Using dimensionless variables, this also reads 



(4) 



The wake force, fwake in equation ([2]), is now modelled using the wake oscillator dynamics proposed in Facchinetti 
et al. [2J| and used for example in Violette et al. [20,l21| or Xu et al. |2j|. The fluctuating lift coefficient q = 2Cl/Clo 
is assumed to satisfy a Van der Pol equation forced by the cylinder acceleration, so that the coupled cylinder and 
wake equations read respectively 



y+ U + -]y + S\j = Mq, 
q + e{q^ -l)q + q = Ay, 



(5a) 
(5b) 



where M = Cio/lB/iTr^S'i^ and A, e are two dimensionless parameters based on experimental data. In all the paper, 
we shall take A — 12, e = 0.3, Co — 2, Clo = 0.8, St = 0.17, ^ — 2.79 and C'mq = 1 as in previous work [20], so that 
M = 0.06 and 7/^ = 0.34 in equations (P. 

Equations ^ are integrated in time using finite differences. Using initial conditions with the cylinder at rest and 
the wake variable subject to a small random perturbation, the system reaches a limit cycle within a few periods. This 
cycle is then analyzed in terms of frequency and time-averaged quantities, such as the efficiency, equation (jj]). 

B. Optimal energy harvesting and lock-in 

In the present case, the harvesting efficiency 77 depends only on the frequency ratio S and the reduced damping ^. 
Figure [2Ka) in fact shows a very strong dependence of the efficiency with the parameters, notably the frequency ratio 
S. The optimal efficiency ry = 0.23 is reached for 6 = 0.89 and £, — 0.20. This is consistent with the value oi rj = 0.22 
obtained experimentally in Bernitsas et al. 18|. The model allows one to recover an expected result which is that the 
efficiency vanishes for both small and large dampings. In the limit of small damping, the amplitude of VIV saturates, 
as is known from the Skop- Griffin plot [7], and the efficiency then varies linearly with ^ and tends to zero. For high 
dampings, the amplitude decreases like 1/^ so that the efficiency, which scales as ^y^, will also vanish. In terms of 
frequency ratio between the solid and fluid frequency 6 = LOs/uf, the optimal is close to one. This corresponds to the 
classical lock-in condition, known to lead to high amplitudes. 





FIG. 2: Effect of the damping ^ and the frequency ratio 5 in the 2D model : (a) Harvesting efficiency using the full non-linear 
model (level step : 0.015), (b) Growth rate derived from the linear model (level step : 0.03) 



To model lock- in, a simplified version of the wake oscillator model was proposed in de Langre [25| . Linearizing 
equations ^, a linear stability analysis of the system was then amenable : lock- in appeared to correspond to the 
highest growth rate of the coupled mode instability involving the cylinder and the wake. In the present case, the 



linear equations read 

^ + - ] y + 6''y = Mq, (6a) 

q — eq + q — Ay. (6b) 

Note that we keep here all the linear terms, including damping, whereas only the dominant terms were used in 
de Langre [25|. Looking for harmonic solutions in time leads to the frequency equation for u 
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+ iuj^ (e - a) + uj^ {ea + AM - 1 ~ S^) + iuj {a - eS^) +6^=0, (7) 

where a = ^ + 7//X. The dynamics of the system will be dominated by the most unstable mode [20jl25|. Its growth 
rate is shown on Figure [2jb) as a function of ^ and S. Figure HJb) displays both the lock-in effect for 6=1 and the 
decrease of the growth rate for large values of ^, two features that were found on the efficiency map (Figure [2k')- In 
terms of efficiency, large values are found to lie in the zone of large linear growth rates. 

This section shows that a simple VIV model with damping allows to capture the main expected effect of the two 
harvesting parameters, the frequency ratio and the intensity of damping/harvesting. The optimal efficiency is reached 
for a clear lock-in condition and a well balanced damping value. 

III. INFINITE CABLE WITH DISTRIBUTED ENERGY HARVESTING 

We consider now an infinite tensioned cable as a three-dimensional generalization of the elastically mounted rigid 
cylinder of the previous section. Again, the harvesting device is modelled by a dissipation introduced in the system, 
here a uniformly distributed damping support (Figure [Dd) . 

The dynamics of the cable and of the wake are modelled as in the previous section, but with the noticeable change 
that the stiffness term hY in equation ([T]) is now replaced by a tension induced stiffness —Qd^Y/dZ^, where Q is the 
uniform tension in the cable and Z the spanwise coordinate [13l |21| . The corresponding dimensionless equations are 

y+U + Ay~y"^ Mq, (8a) 

q + e{q^-l)q + q = Ay, (8b) 

where z = ZiOfy^mt/Q and the derivation with respect to the dimensionless spanwise coordinate z is noted ( ) , 
all other variables and parameters being identical to those of the previous section. Note that we use here a set of 
dimensionless variables that differs from the above cited papers, for the sake of further clarity. These equations are 
now integrated in time and space using finite differences on a spatially periodic domain. Its size is tested to be long 
enough compared with the wavelength of the developing motion for it not to influence the result. Initially, the cable 
is at rest and a small random perturbation is applied to the wake variable q. Only the steady state response is here 
considered. 

The influence of the only harvesting parameter, the damping ^, is now explored. It is found that the system responds 
in the form of propagating quasi- harmonic waves, as previously described in Facchinetti et al. [Ij], Violette et al. [21| 
and illustrated on Figure Efa) . Note that these waves have an amplitude that depends on ^ but their wavelength 
apparently does not. The efficiency is uniform in span and varies as a function of ^ in the classical bell shape (Figure 
|31d). An optimal is found, r/ = 0.22 for ^ = 0.18. These values are very close to those found in the 2D case, but not 
identical. 

To better understand these dynamics, the linearized equations are used, as previously. These read here 

^ + -)y~y" = Mq, (9a) 

q-eq + q = Ay. (9b) 

Looking for a harmonic solution both in space and time, the dispersion relation between the complex frequency cj 
and the real wavenumber k is found as 

u}^ + ioj^ ie-a)+ uj^ {ea + AM - 1 - k^) + iuj {a - ek^) + fc^ = 0, (10) 

where a = ^ -|- 7/a*. This equation is identical to that obtained in the 2D case, equation ([7]), but for the replacement 
of the frequency ratio S = uJs/i-Of by the wavenumber k. We know from the previous section that the highest growth 





FIG. 3: (a) Instantaneous displacement of the cable for different values of the damping for ^ — 0.1 (solid) and ^ = 1 (dashed), 
(b) Comparison of the optimal harvesting efhciency 7720 (^) from the elastically mounted rigid cylinder (dashed) and that of 
the infinite tensioned cable with distributed harvesters, 773D {£,) (solid). 

rate is obtained for S — 1, hence fc = 1 here. This explains two features of the dynamics of the cable. First, the 
numerically-observed wavenumbers are actually within 10~^ of this reference value fc = 1 for all explored values of £, 
(Figure [5^) : the system does always select the linearly dominant wavenumber. This is consistent with the work of 
Violette [20[. As a result, the 3D efficiency curve tj^d (C) (Figure [8)3) is strictly equivalent to the cut of the efficiency 
map of the elastically- mounted rigid cylinder at the level (5=1 (Figure [5^). For each damping ^, the optimal frequency 
ratio on Figure HJa) differs slightly from S = 1. An efficiency curve 772,0 (C) is obtained from the map Figure [^fa) by 
selecting the optimal 6 for each damping. It is plotted on Figure Efb) : the 3D and 2D curves are very close but not 
identical. In particular, the 3D optimal efficiency 77 = 0.22 slightly differs from the 2D optimum 77 = 0.23. 

It is interesting to note that even if the tensioned cable case is not exactly as efficient as the elastically supported 
cylinder one, it has a major advantage in the fact that it adapts its stiffness by adapting its wavelength so that it is 
always at lock-in. A rigid cylinder on the other hand only has a short range of operability in terms of flow velocity, 
which corresponds to tUs = ujf. This will be further discussed in Section FVl 



IV. INFINITE CABLE WITH LOCALIZED ENERGY HARVESTING 



A. Fluid-solid model with discrete energy harvesting 



In practice, most energy extractors are not continuous devices, as assumed in the previous section. We now address 
the case where the corresponding dampers are periodically spaced along the length of the cable (Figure [TJ;). All the 
features of the model are identical to those of the preceding section, except dissipation, which is now governed by the 
the dynamics of each damper. 
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where [ ]_ stands for the jump between values on both sides of a damper and R is the damping coefficient. Here L is 
the distance between two consecutive dampers. We focus on cable motions that are periodic in space over a length 
L, so that the dimensionless equations reduce to 

(12a) 

(12b) 



with the condition 
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where / = Lujf\Jmt/Q is the dimensionless distance between dampers, ^ = R/LmtWf is the equivalent damping 
coefficient per unit length and all other dimensionless variables are chosen as in the other sections. In this case, the 
efficiency reads 
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df 



you^L 



(14) 



Using the variables defined above, the efficiency also reads r/ = 16/i7r'^S'f (^y^)j ^^ in the previous cases. 

Equations ([72]) and ([T3|) are now integrated in time and space, using the same method as in the previous section, 
and the steady state dynamics of the cable is analyzed in terms of efficiency 77. 

B. Optimal energy harvesting 

Two dimensionless parameters, ^ and I, now govern the dynamics of the system. In order to compare with the 
results of the other geometries, we seek, for any given damping ^, the maximum efficiency over all possible length I. 
Figure |3] shows the comparison between the SD-periodic case, the SD-continuous case and the 2D case. It appears 
that (i) the maximum efficiency for the periodic harvesting, 77 = 0.19, is almost as high as for the other cases, (ii) 
the corresponding optimal damping per unit length is much higher than in the previous cases, ^ = 3.65, and (iii) 
the evolution of the efficiency with damping does not follow the classical bell shape. A more detailed analysis of the 
efficiency map is thus needed to understand these features. 
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FIG. 4: EfBciency as a function of the damping coefhcient for the three considered configurations : the 2D case of an elastically 
mounted rigid cylinder (dashed), the infinite tensioned cable with continuous harvesting (solid) and the infinite tensioned cable 
with localized harvesting devices (dash-dot). 



The full map of efficiency is given in Figure [5ja) . The effect of the parameters is complex, and may be analyzed by 
considering three zones. First, at the bottom of the map, or for small distances between harvesting devices (zone A) 
a simple effect of the damping coefficient is recovered : the evolution of the efficiency with the damping coefficient is 
close to the expected bell shape. The global maximum of efficiency also lies in this zone. In Figure [5jb), we also show 
that a large part of the full efficiency curve is explained by the variations of efficiency with the damping in this zone. 
For larger length but low damping, typically ^ < 0.4 (zone B), the efficiency is significantly smaller and only weakly 
depends on the length I. The evolution with ^ is bell-shaped and the maximum over l^ in this zone, corresponds to the 



bottom left part of the full efficiency curve (Figure [51d). A much more complex behaviour is found in the third zone 
(zone C). The efficiency strongly depends on I : tongues of high efficiency are surrounded with negligible harvesting 
regions. These tongues are responsible for the discontinuities in the full efficiency curve as shown on Figure [51(b). 
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FIG. 5; (a) Map of the energy harvesting efficiency from VIV of an infinite cable with discrete harvesting devices as a function 
of the damping density ^ and the distance between two dashpots I and clear identification of the three zones A,B and C of the 
parameters plane and (b) contribution of each zone in the efficiency curve of this configuration of energy harvesting. 



These zones differ essentially by the kinematics of the cable. A particular case is that of zone B, where the cable 
deforms in the shape of propagating waves, of wavelength 27r, Air, 6n, ... depending on the length I. This is illustrated 
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in Figure [6l^a,b) . As the efSciency is low in this zone, we shah not discuss it further. The optimal motion, which lies 
in zone A, is shown Figure ID^c). It is very close to the classical sinusoidal mode 1 of a tensioned cable, but with a 
small displacement of the support, which is responsible for the energy harvesting. All motions in zone A are of this 
type. At the same value of f but for larger distances I, two cases relevant to zone C are displayed in Figure [S{d) and 
(e). For the first one, Figure [6l^d), the efficiency is almost zero and the motion is similar to a mode 2 of a cable, with 
hardly any motion of the support. Conversely, Figure |6l[e) , a high efficiency is achieved with a motion resembling a 
mode 3, with support displacement. 








FIG. 6: Dynamics of the cable in different parts of the efficiency map, the sohd line stands for the instantaneous position of the 
cable, the dashed line for the envelope of its motion and the horizontal dotted line recalls the position of one of the harvesters, 
(a) I = 6.80, C = 0.0193, r/ = 0.05, zone B, (b) I = 13.41, ^ = 0.0193, r] = 0.05, zone B, (c) I = 3.43, C = 3.65, r? = 0.19, zone A 
corresponding with the optimal case, (d) I — 6.80, ^ = 3.65, 77 = 9.10"^, zone C, (e) I = 9.56, ^ = 3.65, rj = 0.07, zone C. 



C. Impact of the mode shapes on eflficiency 

Based on the comparison between the last three cases, we may expect that any motion resembling an odd mode of 
a cable will have a large efficiency, while even modes, which correspond to negligible forcing on the dampers, leads to 
inefficient harvesting. This is illustrated on Figure [71 where the odd mode numbers will naturally exhibit a jump in 
slope, and thereby a force on the dashpot. 
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FIG. 7: Mechanism of excitation of the harvester by stationnary modes, the arrows represent the local tension. Left : odd 
modes, high efficiency. Right : even modes, low efficiency. 



As the mode number seems to be of primary importance in determining the efficiency, the linearized version of 



the equations is used to estimate it, as in previous sections. For each set of parameters {^,1), the most unstable 
mode is analyzed in terms of its wavelength A and a corresponding mode number is defined as n = 21/ X. Note that 
this number may not be exactly an integer, due to the boundary condition : in that case the closest integer is used. 
The regions corresponding to each mode number are reported on the efficiency map (Figure |S]). The regions of high 
efficiency in zone A and the tongues of zone C are all associated with odd mode numbers, as expected. 
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FIG. 8: Correspondence between the mode numbers zones delimited by the solid lines and the efficiency map. 
odd mode numbers regions match very well with the high efficiency zones. 



Note that the 



In zone B, only even modes exist. Odd modes, which are the efficient ones, only appear for high damping, in zone 
A and C. This explains why the optimal efficiency of this configuration is obtained for rather high damping. Indeed, 
damping plays a double role here by transforming the overall mode shape and by controlling the local amplitude at 
the dampers locations. 



V. DISCUSSION 



We may now compare the results of the three configurations analyzed above. The first feature to point out is 
that the optimal efficiencies are comparable, of the order of 0.2. They are achieved with different values of damping 
density ^, but always with a lock-in condition between the wake frequency and a frequency of the solid. These lock-in 
conditions deserve to be discussed further. In the simple 2D case, the lock-in condition corresponds to the coincidence 
between the wake frequency, w/, and the natural frequency of the elastically supported body, Wg. It reads simply 
d = ujs/ujf = 1. In the second case. Section 3, the tensioned cable with continuous damping support has the ability to 
adapt its wavenumber k so that the local frequency of motion synchronizes with the wake frequency, namely fc = 1. 
When another length scale is introduced by separating the dampers with a distance L, the choice of the wavenumber 
is restricted, and lock-in is possible for each mode. Yet, harvesting will be significant only for odd modes, and they 
exist only for high enough damping. Figure El The lock-in condition reads, for these particular modes, / = {2n + 1) tt. 

Although lock-in, as described above, is an essential feature of an efficient harvesting, the optimal conditions differ 
slightly from exact lock-in. We now discuss this in terms of dimensional parameters. The optimal regime for the 2D 
case corresponds to 



D 

27r StU 



0.89, 



C 



rD 

2'!rSt'mtU 



0.18, 



(15) 



10 



and for the tensioned cable with focahzed harvesting devices 



TT f D \ n re ^ (R/L)D , , 

/ \2tt btU J L y rrit ZTTbtmtU 

These show that they are of identical forms, except for the replacement of the rigid cylinder natural frequency ujg in 
P^ by the cable mode frequency 

-c^jJ^- (17) 

L\l nit 

For a given mechanical system, the optimal efficiency is reached only at one flow velocity. Yet, an advantage of the 
tensioned system is that the tension is easier to adapt than a mechanical elastic support stiffness. Noticeably, if the 
tension is caused by the flow, so that 9 ^ [/^, the condition (fT6| over the distance between two dampers, ir/l = 0.92 
i.e. I = 3.43, becomes independent of the velocity. As the velocity varies, the damping ^ may shift, but it has a small 
influence, so that a near optimal efficiency is maintained. 

Using equation (|14p and a value of 77 = 0.18, the power harvested by a cable of total length Ltot — 100 m, diameter 
D — 0.04 m and a 1.5 m/s flow is about 1215 W. If the tension comes from the flow, 6 — pDLtotU^ and fi — 2, 
as in Modarres-Sadeghi et al. [ij|, the dimensional optimal distance between two dampers is about L = 4 m and 
twenty-five dampers are required along the cable to harvest those 1215 W. 

All these results have been obtained under assumptions which shall now be discussed. A first assumption, used 
throughout the paper, is that motion is purely cross- fiow, although VIV dynamics of a cable in flow are known to 



combine both in-line and cross- flow motions 15|. Still, in-line motion amplitudes are generally smaller [28|, and the 
harvesters we consider are purely cross-flow oriented. Assumptions on the cable are also made, by considering it as 
initially straight, and by neglecting non-linear geometrical effects in its dynamics. The wavelength of motion are 
here much larger than the cable diameter, discarding non-linear effects, and tension is such that the natural radius 
of curvature is much larger than the same wavelength. This was the case of cable systems such as in Violette et al. 



2l| or Xu et al. [23|. The fundamental model used throughout the paper is that of a wake oscillator. This approach 
has some limitations, particularly in describing the fine dynamics of the wake, as those computed for instance in 
Bourguet et al. [1^, but it is known to capture well features such as lock- in and generic parameters dependences. As 
all the results given in this paper reduce to simple considerations on those aspects, we expect that they are, at least 
qualitatively, relevant and general enough. 

As a conclusion, we may state that energy harvesting using vortex-induced vibrations seems as feasible and promising 
using tensioned cables as it is using the classical rigid cylinder geometry. It may even have some specific advantages 
such as the ability to harvest energy with large devices, or the ability to permanently adapt the cable dynamics to 
the lock-in condition. 
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